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We study the spectrum and eigenmodes of the QCD Dirac operator in a gauge background given 
by an Instanton Liquid Model (ILM) at temperatures around the chiral phase transition. Gener- 
ically we find the Dirac eigenvectors become more localized as the temperature is increased. At 
the chiral phase transition, both the low lying eigenmodes and the spectrum of the QCD Dirac 
operator undergo a transition to localization similar to the one observed in a disordered conductor. 
This suggests that Anderson localization is the fundamental mechanism driving the chiral phase 
transition. We also find an additional temperature dependent mobility edge (separating delocalized 
from localized eigenstates) in the bulk of the spectrum which moves toward lower eigenvalues as 
the temperature is increased. In both regions, the origin and the bulk, the transition to localization 
exhibits features of a 3D Anderson transition including multifractal eigenstates and spectral prop¬ 
erties that are well described by critical statistics. Similar results are obtained in both the quenched 
and the unquenched case though the critical temperature in the unquenched case is lower. Finally 
we argue that our findings are not in principle restricted to the ILM approximation and may also 
be found in lattice simulations. 

PACS numbers: 72.15.Rn, 71.30.+h, 05.45.Df, 05.40.-a 


I. INTRODUCTION 

By now it is well established that many of the properties of light hadrons are dictated by spontaneous chiral 
symmetry breaking (SySB) due to the nonperturbative structure of the QCD vacuum. However, in the high temper¬ 
ature limit, non perturbative effects are suppressed and chiral symmetry is eventually restored at a certain critical 
temperature. Recent lattice simulations suggest furthermore that both the restoration of chiral symmetry and the 
confinement-deconfinement transition occur at the same temperature Q . A detailed account of the microscopic mech¬ 
anisms leading to SySB and confinement and the eventual restoration and deconfinement transition remains among 
the most formidable challenges in hadronic physics. In recent years the combination of experimental results and novel 
lattice techniques have greatly clarified the nature and properties of these transitions. However we are still far from a 
full microscopic understanding of the chiral phase transition and its relation to the confinement-deconfinement transi¬ 
tion. In this paper we examine how the phenomenon of Anderson localization Q , originally introduced in the field of 
condensed matter physics, may help to understand the mechanism driving the chiral phase transition. Specifically, in 
the context of an Instanton Liquid Model (ILM), we show that near the temperatnre around which the chiral phase 
transition occurs the low lying eigenmodes of the QGD Dirac operator undergo a localization-delocalization transition 
consistent with an Anderson transition (AT) characterized by multifractal eigenstates and critical statistics @]. In 
order to make accessible our findings to a broader audience, we first provide short summaries on the role of instantons 
in the non-perturbative structure of QCD, similarities between the QCD vacuum and a disordered medium, the QCD 
chiral phase transition and finally a brief overview on the properties of a disordered conductor at an AT. In section 
HD we will present our results for the ILM simulations in both the quenched case and the unquenched case for two 
and three flavors. In section ITTTI we discuss the relationship of our results with current lattice simulations followed by 
conclusions in section HU 


A. Instanton Liquid Models and SySB 

Instantons Q are classical solutions of the Yang-Mills equations of motion which minimize the action in Euclidean 
space. They are believed to be the leading semiclassical contribution to the bosonic part of the QCD path integral. 
However the construction of a consistent QCD vacuum based on instantons faces serious technical difficulties. Exact 
multi-instanton solutions are hard to obtain since the Yang-Mills equations of motion for QCD are nonlinear and 
therefore a superposition of single instanton contributions is not itself a solution. Additionally quantum corrections 
may spoil the semiclassical picture implicitly assumed of a QCD vacuum composed of instantons well separated 
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and weakly interacting. These problems have been overcome either by invoking variational principles 0 or by 
phenomenologically fixing certain parameters of the instanton ensemble. The latter case, usually referred to as the 
instanton liquid model (ILM) 0 (for a modern review see j^), yields accurate estimates of vacuum condensates and 
hadronic correlation functions [3 for most light hadrons just by setting the mean distance between instantons to be 
.R « 1 fm (which corresponds to a density iV/y ~ 1 fm“"') and the mean instanton size p w 1/3 fm. Lattice simulations 
have also supported the picture of a QCD vacuum dominated by instantons @. 

One of the main success of the ILM is its convincing explanation of how chiral symmetry is spontaneously broken 
in the QCD vacuum ii The starting point is the fact that the Euclidean QCD Dirac operator has an exact zero 
eigenvalue in the field of an instanton. In an ensemble of (anti-)instantons the overlap of neighboring instantons makes 
the would-be zero modes effectively split around zero. In the limit of large number of (anti-)instantons a continuous 
band spectrum is formed with a spectral density finite at zero. The spectral properties of the low lying modes of the 
Dirac operator are thus controlled by these non-perturbative configurations. It turns out that the order parameter 
associated with spontaneous chiral symmetry breaking (the chiral condensate) is indeed related to the spectral 
density p(e) of the QCD Dirac operator around zero through the Banks-Casher relation |^, 

(^ 1 /;) = - lirn lim , ( 1 ) 

e^O y—i-oo V 

where V is the space-time volume. We are therefore interested in modeling the lowest eigenvalues of the Dirac 
operator. The ILM provides a phenomenological model for these low energy modes of the QCD Dirac operator. In a 
basis of N/2 instantonic zero modes and N/2 anti-instantonic zero modes the matrix elements of the Dirac 
operator take the form (for zero quark mass) 

In the limit of large separations Ry = Zi — Zj between the center of instanton i and anti-instanton j the elements of 
the N/2 X N/2 overlap matrix Tia decay as a power-law, 

~ P^P,/\R^J? (3) 

with Pi and pj their respective sizes and p is the Dirac operator. Many properties of the QCD vacuum can then 
be estimated from the matrix elements of Tja- For example the chiral condensate can be calculated as (ipip) ~ 
— (215MeV)^ 0j@| which is reasonably close to the experimental value. 


B. The QCD vacuum as a disordered medium 

In recent years it has been suggested PIHIIllil that the S^SB in QCD and the phenomenon of conductivity in 
a disordered medium have similar physical origins. Conductivity in a disordered sample is produced by electrons that 
although initially bound to an impurity may get delocalized by orbital overlapping with nearby impurities. Similarly, 
in the QCD vacuum, the zero modes of the Dirac operator though initially bound to an instanton may get delocalized 
due to the strong overlap with other instantons. As a consequence, the chiral condensate becomes nonzero and chiral 
symmetry is broken. 

In the case of electrons the overlap is effective only among nearest neighbors since the electron wavefunction has an 
exponential tail. On the contrary, in the QCD vacuum, quark zero modes can travel long distances in the QCD vacuum. 
Hopping is not restricted to nearest neighbors since the probability for a quark to hop between an instanton and an 
anti-instanton has a power-law tail in the pseudo-particle distance (equation l^l. As we discuss next, eigenvectors of 
quantum systems with this kind of long range hopping are delocalized in four (and higher) dimensions (14LI1 !t|. 

Localization properties of a disordered system are strongly influenced by factors such as the dimensionality of 
the space, the strength of disorder and the hopping range. For short range hopping it is well known that, in three 
and higher dimensions and for a certain range of impurity amount, a system will have a mobility edge separating 
localized from a delocalized eigenstates. Eigenstates around the mobility edge (the Anderson transition region) are 
multifractal and the level statistics (commonly referred to as critical statistics jig) are universal and intermediate 
between Wigner-Dyson (WD) and Poisson statistics (for a review see [13) _ 

An Anderson transition is also observed in disordered systems with power-law hopping HHHIIl provided that 
the exponent of the hopping decay matches the dimension of the space. If the exponent is larger (smaller) than the 
dimension of the space then the eigenstates becomes localized (delocalized). In this case the localization properties 
in the thermodynamic limit depend only on the power-law decay and not on the strength of disorder. Since in QCD 
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at zero temperature the dimensionality of the space (four) is larger thau the hopiug decay power (three, see equatiou 
0 ), oue expects the eigeustates to be delocalized iudepeudeutly of the deusity of iustautous (aualogous to the uumber 
of impurities which determiues the disorder streugth). ludeed iu a receut paper we showed that au effective 
raudom matrix model iucorporatiug the pheuomeuological power-law decay for the matrix elements of the QCD Dirac 
operator in a basis of zero modes describes the spectral correlations of the QCD Dirac operator well beyond the 
Thouless energy Ec Jf^/(|('0^)|L^) (F,r is the pion decay constant and L the length of the system) which sets the 
limit of applicability of standard chiral random matrix theory [T^. However at high temperatures, QCD becomes 
effectively three dimensional, the chiral symmetry is restored and we expect to see localized states appear in the 
spectrum. One of the main aims of this paper is precisely to investigate the relationship between the vanishing of the 
chiral condensate and the localization of the eigenstates of the QCD Dirac operator in a background of instantons. 


C. ILM at finite temperature: chiral phase transition and localization 

In this section we summarize how the ILM is extended to finite temperature and also to what extent reproduces 
the main features of QCD at temperatures around the chiral phase transition. It is clear that for sufficiently high 
temperature (beyond the chiral phase transition) QCD becomes perturbative and instanton configurations should 
be exponentially suppressed. However at temperatures around and below the chiral phase transition the density of 
instantons is still large and it is expected that they play an important role in the description of the phase transition. 
The generalization of the ILM to finite temperature was presented in Below we provide a brief summary of its 
main features and refer to @ for details: 

• The effective coupling constant g{T) decreases (QCD is an asymptotically free theory) as temperature increases. 
The density of instantons p decreases with temperature since p oc exp(—5) with S = 87r^/g^(T). At the 
restoration temperature the density is still roughly half of the zero temperature value. 

• Instanton like solutions (usually called calorons) are still zero modes of the QCD Dirac operator at finite 

temperature. However the decay of the associated fermionic zero mode '^(i?) {R is the 3D instanton- 

anti-instanton distance) is exponential (instead of power-law) in the spatial dimensions and oscillatory in the 
Euclidean time dimension. 

• The probability \Tja\ ~ (a a numerical factor of order unity) for a fermionic zero mode to hop from 

an instanton to an anti-instanton has an exponential tail in the spatial dimensions thus restricting the overlap 
only to nearby instantons. On the contrary, in the time direction the probability is oscillatory. This suggests 
that at finite temperature the isolated fermionic zero modes are exponentially localized in space but delocalized 
in time, thus for localization studies the dimension of the system is three and not four. In this case, since the 
effective hopping is short range, the global localization properties of the eigenstates of QCD Dirac operator are 
expected to depend strongly on the density of instantons as for normal disordered conductors in 3D. 

• The exponentially small tail of the overlap matrix elements together with the oscillatory behavior in Euclidean 
time suggest that in the high temperature limit (higher than the restoration temperature) the instantons and 
anti-instantons are paired in molecules aligned along the time direction. As a consequence the spectral density 
of the QCD Dirac operator develops a gap near zero making, according the Banks-Casher relation, the chiral 
condensate vanish. 

Although certainly appealing, the validity of the ILM at finite temperature as an accurate model for QCD around 
the phase transition should be corroborated by lattice simulations. Recent lattice studies suggest that this is the 
case though further investigations are needed to fully settle this issue. 


D. Anderson transition: description and signatures 

A disordered system with short range hopping in three or more dimensions displays a metal-insulator transition or 
Anderson transition (AT) around the band center for a critical amount of disorder. By critical amount of disorder we 
mean if the disorder is increased all the eigenstates become exponentially localized. For a disorder strength below the 
critical one, the system has a mobility edge at a certain energy which separates localized from delocalized states. Its 
position moves away from the band center as the disorder is decreased. Delocalized eigenstates, typical of a metal, 
are extended through the sample and the level statistics agree with the random matrix prediction for the appropriate 
symmetry. In this paper we will focus on QCD with three colors which corresponds to random matrix predictions from 
the Gaussian Unitary Ensemble (GUE) |^. The level statistics obtained from random matrix theory are also referred 
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to as Wigner-Dyson statistics (WD). On the contrary, localized eigenstates, typical of an insulator, are characterized 
by an exponential decay \'4>{t)\ ~ with the localization l eng th ^ smaller than the system size. The spectral 

correlations in this case are well described by Poisson statistics corresponding to uncorrelated eigenvalues. 

Around the mobility edge the system is referred to as being at the AT (or at the metal-insulator transition). 
Signatures of this transition are found in both the spectrum and the eigenfunctions. The eigenstates are multifractal, 
namely, the wavefunction moments Pq = J d‘^r\’tp{r)\'^'^ oc present anomalous scaling with respect to the 

sample size L, where Dq is a set of different exponents describing the transition |^. The second moment P 2 
is usually referred to as the inverse participation ratio (IPR). Despite the fact that multifractal eigenstates occupy 
a vanishing volume, they strongly overlap each other. For instance, the density-density correlation function of two 
eigenstates with energies and decay as a power-law, 

J d^r\^nir)\^\^ljmir)\^ ~ \Er, - (4) 

(for a review see [iTil^ b 

Level statistics at the AT are universal and intermediate between WD and Poisson statistics. Typical features 
include: 

1. A scale invariant spectrum such that any spectral correlator utilized to describe the spectral properties does 
not depend on the system size. We recall that, according to the one parameter scaling theory of localization, 
the AT occurs when the disorder strength is such that the conductance is size independent (the beta function 
has a zero). 

2. Level repulsion such that the nearest neighbor level spacing distribution, P{s), (which gives the probability of 
having two eigenvalues at a distance s) vanishes as s —> 0, thus adjacent levels still repel each other as in metallic 
samples. However P{s) ~ for s 3> 1, similar to the result for an insulator P{s) = though the constant 
A is always larger than the unity (~ 1.7 at the 3D Anderson transition). 

3. Long rang e correlators such as the number variance S^(^) = {{Ne — {Ni))^) ^ for ^ 1 are asymptotically 

linear |^, as for an insulator S^(^) £, but with a slope X < 1 (y 0.27± .02 for a 3D Anderson model). Here 

Ni is the number of eigenvalues in an interval of length i. Likewise the spectral rigidity (a smoothed version of 
the number variance) given by = -A (£3 _ _|_ x^)Y?{x)dx is also linear with slope Aa).^) ~ (y/lS)!*. 

The above discussion, with small modifications, holds independently of the details of the microscopic disorder 
distribution or whether the system has any additional symmetries. Our case of interest, the QCD Dirac operator, 
has an additional chiral symmetry which makes the matrix representing the operator have a block structure given in 
equation ©• Obviously the additional chiral symmetry does not affect the properties in the bulk of the spectrum 
but certainly may have an impact on the eigenvalues close to the origin. It turns out that in this region |lM 1^ 1^ 
the properties of the eigenstates depend strongly on the dimension of the space and the disorder distribution. For 
instance for a 3D disordered system with chiral symmetry and short range (nearest neighbor) disorder |2(il 1^ the 
eigenvectors are power-law localized with and exponent depending on the density of impurities. Typically, in lower 
dimensions (ID,2D) eigenvectors close to the origin are more delocalized than those at the bulk. The reasons for 
such anomalous behavior have been traced back to the absence of weak-localization corrections in the chiral case. 
By contrast in three dimensions, it has been claimed that eigenstates close to the origin are more localized than 
those in the bulk for a broad range of disorder strengths. However there is no consensus on whether this finding is 
universal of it is just a peculiarity of the model utilized in |2 (t|. In any case it is expected that, for a certain disorder 
strength, an AT close to the origin may coexist with the one above reported at the bulk. Indeed, this is also what we 
have found in the ILM at finite temperature (see next section). Clearly further investigation is still needed to provide 
a detailed account of the AT close to the origin in disordered systems with chiral symmetry. 

In summary, at zero temperature certain features of the QCD vacuum resemble those of a disordered medium where 
the role of impurities is played by an instanton liquid. The SySB is a consequence of the delocalized nature of the 
zero mode eigenvectors of the QCD Dirac operator caused by the long range nature of the overlap matrix. A natural 
question to ask is whether these analogies can be extended to the case of finite temperature. The aim to this paper is 
to answer affirmatively this question. Specifically, in the context of an Instanton Liquid Model, we show that around 
the temperature that the chiral phase transition occurs the low lying eigenmodes of the QCD Dirac operator undergo 
an Anderson transition, characterized by multifractal eigenstates and critical statistics |3 • 
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II. ANDERSON LOCALIZATION IN THE ILM AT FINITE TEMPERATURE 

In this section we show that both the spectrum and the eigenvectors of the QCD Dirac operator in an ILM undergo 
a transition to localization at a certain temperature. Both the quenched and the unquenched cases are considered 
here. In the latter case we see that the transition to localization occurs around the same temperature as the chiral 
phase transition thus suggesting that both phenomena are deeply related. This is the main result of this paper. First 
we provide a short overview on how the numerical simulations were performed. 


A. Technical details of the numerical simulation 

The ILM partition function for Nf quark flavors with masses is given by 

r 

Zinst = Dn n (5) 

fc=i 

where the integral is over the positions, sizes and orientations of the instantons and S'ym is the Yang-Mills action. The 
fermion determinant is evaluated in the space of the fermionic zero modes of the instantons. For further discussion 
of this partition function we refer to 0. We just mention that we use the standard phenomenological value of the 
instanton density N/V = lfm~^. For the sake of simplicity we have kept this density fixed as the temperature 
is increased. We justify this approximation based on the fact that even for temperatures close to the chiral phase 
transition the instanton density is still sizable (around 0.6 fm~'^ [^). Indeed in previous simulations the drop in the 
density was ruled out as the physical mechanism leading to the chiral phase transition [^. 

All units in the ILM are typically given in terms of the QCD scale parameter A which we have simply set to 200 
MeV. This choice is sufficient for our purposes since we are not concerned with making quantitative predictions about 
the position of the QCD chiral phase transition, but instead are interested in more general features of the transition 
that we don’t expect to be very sensitive to a fine tuning of the parameters of the ILM. We therefore use the ILM 
as a qualitative model for QCD at temperatures around the chiral phase transition. We also stress that it provides 
a reasonable description of S^SB and many hadronic correlation functions both at zero and at finite temperature. 
Furthermore the ILM has modest computational requirements and allows us to go to fairly large volumes and get 
good statistics. 

In this paper we present results for both the quenched Nf = 0 and the unquenched Nf = 2, 3 cases {Nf = 2 with 
rUu = md — 0 and at Nf = 3 with m„ = rud — 20 MeV and iris = 140 MeV). The partition function, equation 
is evaluated using a standard Metropolis algorithm. In the quenched case we performed between 2000 to 10000 
measurement sweeps for each set of parameters after allowing around 1000 sweeps for thermalization. Results are 
presented for ensembles of up to V = 5000 instanton and anti-instantons. In the unquenched case we were only able 
to investigate ensembles of up to V = 700 instantons and anti-instantons. We also discarded the first 1000 sweeps in 
each simulation and did between 1000 to 5000 sweeps per ensemble. Thermodynamic properties of the model have 
already been addressed in and will not be discussed here. We will focus mainly on observables such as level statistics 
and eigenvector sc aling properties of the QCD Dirac operator which are especially well suited for studying Anderson 
localization effects . All the spectral correlators are calculated from the unfolded spectrum. This procedure scales 
the eigenvalues so that the spectral density on a spectral window comprising several level spacings is unity but it does 
not remove the small fluctuations about the mean density that provide fundamental information about the system. 

Finally we recall that in the ILM the effective disorder parameter is the temperature since we have fixed the density 
of instantons. One of the most challenging tasks of the numerical calculation is to accurately locate the temperature 
at which the AT occurs. In all cases in the paper this temperature has been estimated by the finite size scaling 
method introduced in in the context of disordered systems. In essence this method consists of computing a 
spectral correlator, such as the level spacing distribution P{s) or the number variance for different sizes and 

then finding the temperature at which it becomes size independent. 


B. Anderson localization in the ILM: The quenched case 

In this section we present results for the quenched ILM at nonzero temperature. Without dynamical fermions there 
is not technically a chiral phase transition since there is no chiral symmetry. However one can in principle study 
the “quenched” quark condensate in a purely gluonic background through the Banks-Casher relation 0. Lattice 
simulations with staggered fermions Q have also found that the vanishing of the “quenched” quark condensate 
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FIG. 1: Inverse participation ratio (IPR) times N versus the 
average eigenvalue for the quenched ILM with N = 2000 in- 
stantons at different temperatures. Low lying eigenmodes are 
more localized than those in the bulk for all tested tempera¬ 
tures. 


FIG. 2: Level spacing distribution P{s) for the quenched ILM 
with N = 5000 instantons at T = 200 MeV for different spec¬ 
tral windows. The legend indicates the eigenvalue numbers 
with 0 being the smallest one. The spectral correlations of 
the lowest and highest modes are close to Poisson statistics. 
By contrast the spectral correlations in the central part of 
the spectrum are well described by RMT (GUE result). 


occurs at the same temperature as the deconfinement transition does. Since the ILM does not have a confinement- 
deconfinement transition 0 we can only speak about the chiral phase transition in the ILM. However, in the quenched 
ILM, the spectral density seems to diverge close to the origin (similarly to recent lattice results with overlap fermions 
1^1 thus suggesting a likewise divergent “quenched” quark condensate. In view of these facts the transition to 
localization (see below) we have observed in the quenched QCD Dirac operator at finite temperature both close to 
origin and in the bulk of the spectrum can’t be linked directly to either of these phenomena. We still consider the 
quenched ILM results to be of interest as we have found a very clear example of a mobility edge with features strikingly 
similar to those of a 3D disordered system at the AT. This finding reinforces the idea [I^ that certain aspects of the 
QCD vacuum resemble those of a disordered system. 

We start by summarizing our main findings: 

1. The eigenvectors of the QCD Dirac operator at finite temperature with the lowest eigenvalues (the origin) are 
much more localized than the ones in the middle of the spectrum (the bulk). Signs of stronger localization close 
to the origin were observed not only in the eigenvector moments (figure^ but also in the level statistics (figure 
I2J). This difference between the origin and bulk has also been reported in 3D disordered systems with chiral 
symmetry |2 (t). 

2. In the bulk we hnd a mobility edge separating localized from extended states which moves toward the high end 
of the spectrum as the temperature is decreased. The level statistics of eigenstates around the mobility edge 
are well described by critical statistics (see figures 13 and 0)1 typical of a 3D disordered conductor at the AT. 
Outside the critical region the level statistics are similar to those of a metal (Wigner-Dyson statistics) in the 
region below the mobility edge. By contrast, in the region above the mobility edge, it resembles those of an 
insulator (Poisson statistics). 

3. Close to the origin we observe a transition to localization for temperatures in the range T ~ 100 — 140 MeV 
(see figure 13 . 


1. The bulk 


Using the finite scaling method introduced in we found a mobility edge in the bulk of the spectrum in the range 
T ~ 150 - 250 MeV. As the temperature decreases its location moves to the end of the spectrum. For T < 170 MeV 
the results are less reliable since the mobility edge is located almost at the end of the spectrum where truncation 
effects due to the finite size of the ILM are larger. Following the literature in disordered systems we have investigated 
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FIG. 3: Level spacing distribution P{s) in the bulk of the 
spectrum for the quenched ILM at T = 200 MeV and different 
numbers of instantons. The inset shows the tail of P{s). The 
solid line corresponds to the best fit P(s) ^ 



FIG. 4: Spectral rigidity A 3 (^) in the bulk of the spectrum for 
the quenched ILM at T — 200 MeV and different numbers of 
instantons. The result is almost size independent and agrees 
well with the prediction of critical statistics (solid line). 


the temperature T ^ 200 MeV such that the mobility edge is located around the center of the spectrum. We initiate 
our analysis by looking at the level statistics. We recall the spectral correlators were computed from the unfolded 
spectrum in the region around the mobility edge. The main findings are summarized as follows: 

1. The spectrum is scale invariant to high degree. As observed in figured the level spacing distribution P(s) does 
not depend much on the system size for volumes ranging from N = 500 to N = 5000. The inset plot shows that 
even the tail of P{s) which contains information about correlations at larger scales is also scale invariant. The 
study of long range correlators further confirms this behavior. Likewise, in figure^] it is shown that the spectral 
rigidity varies little with the system size in a window of more than 200 eigenvalues around the mobility edge. 

2. Level repulsion, P{s) ^ 0 as s ^ 0, typical of a metal is still present (figure^). However P{s) ~ g-i-86s 
exponential tail as for an insulator. Both features also appear in the case of a 3D disordered system at the AT. 

3. The spectral rigidity (figure 0]) is asymptotically linear with a slope y/ 15 ^ (0.29 ± 0.02)/15 in fair agreement 
with the result of a 3D disordered system at the AT of y ~ 0.27/15 

4. The level statistics of the ILM in the critical region are accurately described by generalized random matrix 

models capable of reproducing critical statistics (see figure QJ . The explicit analytical result for the 

chiral case can be found in |32| for a random banded model with a power-law decay and in |3lL 1^ in the 
context of the Calogero-Sutherland model at finite temperature. We recall that in these models there exits 
a free parameter that must be fixed by, for instance, matching the slope of the spectral rigidity y. In our case, 
as shown in figure^ the exact form of the spectral rigidity in the bulk follows closely the prediction of critical 
statistics, equation (31) in with h, a free parameter, set to a value of 0.62. 

5. The level statistics outside the AT window are given by Wigner-Dyson statistics (typical of a metal) in the region 
below the mobility edge. By contrast it is close to Poisson statistics typical of an insulator in the region above 
the mobility edge. This is corroborated in figure [21 where we observe a transition from Poisson to Wigner-Dyson 
statistics in P{s) at T = 200 depending the region of the spectrum where P{s) is computed. 

The above analysis provides compelling evidence that for a certain range of temperatures the level statistics in the 
bulk of the QCD Dirac operator in the ILM are consistent with those of a 3D disordered system at the AT. 

We now turn to eigenvector properties. As mentioned previously, a signature of the AT is the multifractality of 
the eigenstates. Multifractality is usually detected by the anomalous scaling of the moments of the eigenfunction Pq 
(defined in section ITDli with the system size L oc . We have computed the fractal dimension £>2 by fitting P 2 
for different system sizes N. A fit to the form P 2 = aN^ yields b ~ 0.50(5) and consequently D 2 ~ 1.5(1), consistent 
with the D 2 ~ 1-4(2) found for a disordered system at the 3D AT. In order not to consider eigenstates outside the 
critical region we have taken only 10% of the eigenvectors around the center of the spectrum. 















FIG. 5; Level spacing distribution P{s) close to the origin of 
the spectrum for the quenched ILM with N = 2000 instantons 
for different temperatures. Agreement with the prediction of 
RMT (GUE result) is observed for relatively low tempera¬ 
tures T < 100 MeV. A transition between RMT (GUE) and 
Poisson statistics is observed in the range T ~ 100—140 MeV. 
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FIG. 6: Number variance E^(i!) at the origin for the quenched 
ILM with N = 1600 instantons for different temperatures. 
Agreement with the RMT (chiral GUE) result is observed for 
relatively low temperatures T < 80 MeV. The temperature 
range of the transition to Poisson agrees with that obtained 
from the level spacing distribution. 


To recapitulate, the numerical analysis shows with great accuracy that at T = 200 MeV the QCD Dirac operator 
in a background field given by the quenched ILM undergoes and metal-insulator transition in the central part of the 
spectrum. The eigenstates around the transition region are fractal and the level statistics are well described by critical 
statistics, namely, they have all the properties expected in an AT: scale invariance, level repulsion and sub-Poisson 
spectral rigidity. Furthermore the numerical value of parameters such as the fractal dimension D 2 and the slope of 
the number variance y are very close to the results for a standard 3D disordered system at the AT. 


2. The origin 

This region is of special interest since the smallest eigenmodes are responsible for SySB. However as mentioned 
earlier, in the quenched case there is not, in a strict sense, chiral symmetry so the localization properties in this 
region cannot be related neither to SySB nor its restoration at finite temperature. Despite these limitations we have 
also investigated both the spectrum and the eigenmodes in this region. Our main motivation is to test whether the 
similarities between the QCD vacuum (as given by the ILM at a certain temperature) and a 3D disordered system 
at the AT can also be extended to region close to the origin thus reinforcing the idea put forward in that both 
models belong to the same universality class. 

As in the bulk, the level statistics are temperature dependent. For low temperatures T < 80 MeV (see figures 0 
and the eigenmodes are delocalized and the level statistics are very close to the RMT prediction as given by the 
GUE. In the region T ~ 100 — 140 MeV we observe a transition to localization. By contrast at higher temperatures 
the spectrum is very nearly Poisson like as for an insulator. 

We remark that the numerical results close to the origin are less conclusive than those in the bulk. On the one 
hand the analysis is complicated by the accumulation of very small eigenvalues present in the quenched ILM. On 
the other hand the statistics are worse since only the smallest eigenvalues are affected by the chiral symmetry of 
the model. Despite these technical problems it seems clear that the QCD Dirac operator undergoes a transition to 
localization for T ~ 120 MeV. As shown in figure Q] eigenvectors for the smallest eigenvalues, due to the additional 
chiral symmetry, are more localized than those in the bulk. This has already been observed in lattice simulations 
[ 23 , and 3D disordered systems with chiral symmetry . We do not have a clear understanding of the reason for 
such behavior. Indeed ID and 2D chiral disordered models usually show the opposite behavior, namely, states in the 
bulk are localized for any amount of disorder and only the states close to the origin (the ones affected by the chiral 
symmetry) are truly delocalized. 
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FIG. 7: IPR times the total volume (P = in units of fm) 
versus the average eigenvalue for the ILM with 2 massless 
quarks at different temperatures and spatial volumes (L®, in 
units of fm^). 



FIG. 8: Chiral condensate ('i/’V’) a-nd the IPR of the lowest 
eigenmode versus temperature for the ILM with 2 massless 
quarks at different spatial volumes (L®, in units of fm®). Re¬ 
markably the drop in the condensate occurs in the same range 
of temperatures as the transition to localization of the lowest 
eigenmodes of the QCD Dirac operator. 


C. Anderson transition and chiral phase transition: The unquenched case 

In this section we investigate the eigenvalue and eigenvector statistics of the QCD Dirac operator with two massless 
quark flavors and for three flavors with = rud = 20 MeV and TOs = 140 MeV. In the massless case a true chiral phase 
transition is supposed to occur at a certain (chiral restoration) temperature. Our main aim is to investigate whether 
the chiral phase transition and the localization transition in the QCD Dirac operator occur at the same temperature. 
Below we show that this is indeed the case thus suggesting that the phenomenon of Anderson localization may be the 
fundamental mechanism driving the chiral phase transition. 

We mention that in the massive case there is not a chiral phase transition but instead a crossover as is believed 
to occur in full QCD. The quark masses have been chosen larger than the experimental values in order to observe 
this crossover in the ILM. Additionally such large mass values are closer to the ones typically used in current lattice 
calculations. 

The behavior of the ILM at finite temperature for different numbers of flavors has been studied previously [^ . 
The case of just one flavor Nj = 1 is special since the chiral symmetry is explicitly broken by the U{1) anomaly, 
a condensate is formed but is not related to the spontaneous breaking of any symmetry. This condensate decreases 
smoothly with the temperature and there is no sign of any phase transition [^. By contrast for Nf > 5 the chiral 
symmetry is no longer broken (in the chiral limit) at zero temperature. This is due to the fact that the fermionic 
determinant in the partition function suppresses configurations with low-lying eigenvalues. The rate of suppression 
increases with the number of flavors as . It is believed that for Nf ~ 5 the spectral density of the Dirac operator 
develops a gap close to the origin and, by the Banks-Casher relation, the chiral symmetry is not broken even at zero 
temperature. 

We start our analysis with the eigenmode properties of the QCD Dirac operator with two massless flavors. In figure 
[7|we show the IPR times the total volume V (or number of instantons N) versus the average eigenvalue for different 
temperatures and spatial volumes. Completely extended states will have a constant IPRxM independently of the 
system size, while for completely localized states IPRxR will be proportional to the volume. As in the quenched 
case eigenstates close to the origin are more localized than those in the bulk but the difference between the two 
regions is smaller than in the quenched case. This is not surprising since the fermionic determinant suppress low-lying 
eigenvalues thus weakening the chiral symmetry and consequently the difference between the origin and bulk. At the 
lowest temperature (T = 100 MeV) we clearly see that the scaled IPR is independent of the volume and consequently 
the eigenmodes are delocalized. The opposite occurs for T = 150 MeV where scaling is observed in all eigenmodes, 
especially the ones close to the origin, indicating that the eigenmodes are localized. An intermediate situation happens 
at T ~ 120 MeV where we see some volume dependence near the origin and at the high end of the spectrum while the 
bulk states remain fairly constant. This indicates that a mobility edge separating localized from delocalized states is 
forming at the origin. A natural question to ask at this point is whether this transition to localization occurs in the 
same range of temperatures as the chiral phase transition. 
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FIG. 9: Number variance at the origin for the ILM with two massless quark flavors at different temperatures and spatial 

volumes (L^, in units of fm^). The scaling with volume is very small for T ~ 120 MeV which indicates that the critical 
temperature associated with the AT is very close. 


In figurelHlthe chiral condensate versus temperature is shown for a range of system sizes. We see that the condensate 
decreases as the temperature is increased in a manner consistent with a second order phase transition. Although we 
would need to go to much larger volumes to determine the transition temperature more accurately there is a clear 
bend in the condensate around a temperature of Tc ~ 120 MeV which indicates that the chiral phase transition should 
occur around this region. Remarkably, as shown in figures 0 and 0 this is also around the temperature at which 
the IPR of the lowest eigenmode begins to rise signaling a transition to localization. This finding suggests that both 
transitions could be intimately related. Indeed it shows that the phenomenon of Anderson localization may be the 
fundamental mechanism leading to the chiral phase transition in QCD. In order to further confirm that this novel 
relation between QCD and the theory of disordered system is not a peculiarity of the ILM utilized in this paper we 
are currently carrying out similar numerical calculations in the framework of lattice gauge theory . We stress that 
the relation between an Anderson transition and the chiral phase transition is the main result of this paper. 

We now turn to the level statistics. In order to locate the localization transition temperature more precisely we look 
at the number variance near the origin. In figure 0 we plot this for different temperatures and spatial volumes. 

At T = 100 MeV, the number variance decreases as the system size increase thus indicating that the states are still 
extended at this temperature. At T = 120 MeV there is still a big difference between the smallest volume = 63 
fm^) and the larger ones, but the two larger volumes fall directly on top of each other. We believe that this is very 
near the localization transition and that the smaller volume is just too small to reliably see the proper scaling. For 
higher temperatures T ~ 150 MeV the number variance moves rapidly towards the prediction of Poisson statistics 
11 ^(1') = £. typical of an insulator. 

We observe a similar picture for short range spectral correlators such as P{s). For a fixed volume the spectral 
correlations diminish as we increase the temperature. Thus the level spacing distribution moves (see figure irnil from 
the RMT result (GUE) typical of a metal to Poisson statistics typical of an uncorrelated spectrum. By fixing the 
temperature at 120 MeV and looking at the spatial volume dependence i figure [TT]| of P(s), we observe very little 
difference for the range of sizes considered in agreement with the results for the number variance. We also see that 
there is still level repulsion but the tail of P{s) is exponential as for an insulator. All these features are typical 
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FIG. 10: Level spacing distribution P{s) at the origin for the 
ILM with 2 massless quarks with spatial volume L® = 189 fm^ 
for different temperatures. A transition from the RMT pre¬ 
diction (GUE) typical of a metal to the Poisson result typical 
of an insulator is observed as the temperature is increased. 



FIG. 11: Level spacing distribution P{s) at the origin for the 
ILM with 2 massless quarks at the critical temperature Tc = 
120 MeV for different spatial volumes. The scale invariance 
indicates the system is undergoing an AT. 


signatures of an AT. 

The AT is also characterized by the multifractality of the eigenstates. We investigate this feature by looking at the 
scaling of Pq ~ for low eigenvalues at T = 120 MeV for system sizes ranging from = 63-252 fm^. The 

resulting fractal dimensions Dq are shown in tableOlfor an eigenvalue of (A) = 0.076 which is the smallest eigenvalue 
for the = 63 data. The IPRs for the remaining sizes were interpolated to this eigenvalue (see figure 0 . We 
clearly see Dq change with q thus showing that the eigenstates are truly multifractal. It is worth mentioning that the 
numerical value of Dq is very sensitive to small temperature changes. Thus at T = 111 MeV we find that D 2 = 3.07(5) 
like for a metal {D 2 = 3) but at T = 150 MeV, D 2 = 0.96(4) which is closer to the value for an insulator {D 2 = 0). 

Just in passing we mention that a mobility edge and the corresponding AT has also been found in the bulk of the 
spectrum. As in the quenched case, the temperature at which it appears in the central part of the spectrum is higher 
than the one close to the origin. This is consistent with the picture that, at a fixed temperature, eigenvectors close to 
the origin are more localized than those in the bulk of the spectrum. Indeed similar results have also been obtained 
in lattice calculations 1^ . The eigenvalue and eigenvector statistics are also similar to the results shown for the 
quenched case at the AT so we will not present them here. 

Finally we discuss the case of the Nf = 3 ILM with masses m„ = nid = 20 MeV and rus = 140 MeV. As a general 
comment we observe strong similarities with the quenched and Nf = 2 cases studied previously. Thus, as the plot of 
the level spacing distribution P{s) in figure IT^ shows, the lowest eigenmodes move from extended (metal) to localized 
(insulator) as the temperature is increased. 

More importantly, we have also found that the the localization transition of the low-lying eigenmodes of the QCD 
Dirac operator occurs around the same temperature Tc ^ 125 MeV as the crossover to the chirally symmetric phase 
(which is also not far from the value for two massless flavors). However we remark that since the lightest masses 
were still fairly large the eigenvalue density increases at the origin like in the quenched case and in lattice simulations 
with overlap fermions . Even though the chiral condensate only signals a crossover we still see evidence of critical 
statistics in the critical region Tc ~ 125 MeV. Thus, as shown in hgure^l the level spacing at the origin is scale 
invariant and exhibits level repulsion at small distances while having an asymptotically exponential tail. 


q 

2 3 4 5 

Dq 

2.67(4) 2.34(5) 2.08(6) 1.87(7) 


TABLE I: Fractal dimensions (Dq) for the ILM with 2 massless flavors at T = 120 MeV. 
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FIG. 12: Level spacing distribution P{s) at the origin for 
the ILM with N = 300 instantons and 3 massive flavors at 
different temperatures. As in previous cases a transition to 
localization is observed as the temperature is increased. 



FIG. 13: Level spacing distribution P{s) at the origin for the 
ILM with 3 massive flavors at T = 125 MeV for different 
spatial volumes. The spectrum is nearly scale invariant as 
expected at an AT. 


III. DISCUSSION AND RELATION TO LATTICE RESULTS 

In this section we investigate to what extent the relation between Anderson transition and chiral phase transition 
reported in previous sections is a particularity of the ILM or could be a real feature of the strong interactions. 
Obviously a conclusive answer to this question can only be given after a full analysis of the QCD Dirac operator 
around the restoration temperature is carried out in a direct simulation of QCD on a lattice. 

However it is instructive to clearly establish the reasons for the appearance of the AT in the ILM. As is known, the 
fermionic zero modes in the field of an instanton at nonzero temperature have an exponential tail e“''^ in the spatial 
directions and are oscillatory in the time direction |^. This suggests that the overlap among different zero modes 
is essentially restricted to nearest neighbors in the spatial directions. However, in the time direction, different zero 
modes strongly overlap due to the oscillatory character of the eigenmodes. This situation strongly resembles a 4D 
disordered conductor in the tight binding approximation (only nearest neighbor hopping) with one dimension (time) 
much smaller than the rest so the system can be considered effectively three dimensional. It is well established that 
such a system may undergo an AT depending on the disorder strength. In our case the disorder role is played by 
temperature since the low-lying eigenvectors of the QCD Dirac operator decay as ~ From this discussion it is 

natural to find an Anderson transition in the ILM for a particular value of the temperature. The principal ingredient 
to reach the AT is an exponential decay of the eigenmodes explicitly depending on the temperature together with the 
possibility to tune the effective range of the exponential in order to reach the transition region. Any theory with these 
features very likely will undergo an AT for some value of the parameters even if the objects responsible for localization 
are not the classical instantons. In fact this is not the only scenario in which an AT could arise. As mentioned in the 
introduction, an AT is also expected in the case that at a certain temperature the dominant gauge configurations are 
such that the resulting low lying eigenmodes are power-law localized |^/’(r)| ~ l/r“ with a ~ 3. 

We also remark that the AT in the ILM cannot be induced by instanton-anti-instanton molecules. As the tempera¬ 
ture is increased overlapping among different instantons becomes more rare and eventually, in the the high temperature 
limit, the ILM is composed of independent instanton-anti-instanton molecules. However such situation does not occur 
at temperatures around the chiral phase transition due to the multifractal character of the eigenstates. In the molec¬ 
ular phase the zero modes are attached to a unique instanton-anti-instanton pair, the wavefunction will be localized 
in between the two pseudoparticles and consequently the eigenstates cannot be truly multifractal. At most Dq may 
be a constant not depending on the moment q. Another indication that the molecular phase has not been yet reached 
at the restoration temperature is the fact that the chiral condensate decreases smoothly over a quite broad range of 
temperatures. If the molecular phase is dominant a gap appears in the spectral density and the condensate vanishes. 

We now discuss the relation of our findings with previous lattice results. Lattice investigations of the QCD Dirac 
operator at finite temperature have attracted considerable attention in recent years. Results close to the restoration 
temperature qualitatively agree with the ones reported in this paper. It has been observed that eigenvectors 

close to the origin are more localized than those of the bulk (see figures ^ and ITjl. The eigenvalue density of the Dirac 
operator in the quenched ILM more closely resembles that of lattice simulations of overlap fermions than those 
of staggered fermions [^. This is not surprising since the ILM is based on the low energy topological modes that 
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overlap fermions are supposed to treat accurately while staggered fermions (at large lattice spacing) do not. For two 
massless flavors the ILM does not have a large eigenvalue density near the origin and we can clearly see the decrease 
in the condensate as the localization transition occurs. We cannot make a direct comparison to lattice simulations 
here due to the large cost of simulating massless quarks on the lattice. We also note that a mobility edge has been 
observed in the low energy modes of Wilson fermions |38| however they do not have a chiral symmetry and we cannot 
make a direct comparison. In order to further clarify this issue, an updated detailed lattice analysis of level statistics 
and eigenvectors in the framework of lattice gauge theory would be highly desirable. 

Finally we mention that the true nature of the topological objects responsible for SySB and confinement has not 
been fully settled in lattice QCD. Although there is some support for instanton-like objects there are other studies 
that suggest a different picture. Recently it has been argued that the topological charge may not be localized in 
four dimensional regions of space-time but rather spread out over a lower dimensional manifold |39| . The resulting 
low-lying eigenmodes of the QCD Dirac operator are still delocalized and fractal as we observe in the ILM. This is 
also not far from the results of the ILM in 1I1II3 where it was found that the low-lying eigenmodes are delocalized 
in the sense that in the thermodynamic limit the level statistics are described by random matrix theory. It would be 
interesting to compare these predictions with results from overlap fermions at zero temperature. 


IV. CONCLUSIONS 

We have studied the localization properties of the eigenvalues and eigenvectors of the QCD Dirac operator in a 
instanton liquid model at finite temperature. Due to the chiral symmetry of the model, eigenstates close to the 
origin are more localized than those near the center for any temperature. Near the origin we found a clear transition 
from delocalized to localized states as the temperature is increased. There is a also mobility edge in the bulk of the 
spectrum that moves toward the region of lower energies as the temperature is increased. In both cases, around this 
mobility edge, the eigenvectors are multifractal and spectral correlations are well described by critical statistics as 
for a disordered system undergoing an AT. Remarkably both the transition to localization close to the origin and the 
chiral phase transition signaling the restoration of the chiral symmetry occur at roughly the same temperature. This 
suggests that the phenomenon of Anderson localization may be the fundamental mechanism driving the chiral phase 
transition. 
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